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Abstract. A review of recent progress in numerical studies of the Anderson transition in 
three dimensional systems is presented. From high precision calculations the critical exponent 
V for the divergence of the localization length is estimated to be = 1.57 ± 0.02 for the 
orthogonal universality class, which is clearly distinguished from v = 1.43 ± 0.03 for the 
unitary universality class. The boundary condition dependences of some quantities at the 
Anderson transition are also discussed. 
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1 Introduction 

It is more than four decades since Anderson pointed out the existence of a disorder 
induced metal-insulator transition (the Anderson transition) [|] . After the discovery of 
the scaling behavior of coductance 1^, || , it was realised that the transition is a second 
order phase transition. This naturally leads us to expect a degree of universality, i.e., 
the critical behavior does not depend on the details of the model but depends only on 
the basic symmetries of the system such as time reversal symmetry (TRS). 

Numerical calculations have played a very important role in the investigation of 
random systems, for example percolating systems or spin glasses, where a quantitative 
understanding cannot be obtained with analytic methods. This is also true of the 
Anderson transition. In this paper, we review recent progress in numerical studies of 
the Anderson transition. Special emphasis is put on the universality of the exponent 
V for the divergence of the typical length scales. We discuss the universal behavior of 
the statistics of energy levels and conductance at the transition. We finally comment 
on the bomidary condition dependence of certain quantities such as level statistics, 
the conductance distribution function, and localization length. 
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2 Critical exponents 

The Anderson transition is characterized by a vanishing conductivity, a diverging di- 
electric constant, and a divergence of the length scales as we change the parameter x 
such as Fermi energy (electron density), impurity concentration, and pressure. Near 
the critical value Xc, the vanishing of conductivity a is described by the critical expo- 
nent s 

C^~(x-Xc)^ (1) 
and the divergence of the dielectric constant e 

(2) 



{Xc - xY' ' 

Here we suppose that x > {<)xc is the metallic (insulating) regime. In the metallic 
regime, the characteristic length scale ^ is the correlation length, while in the insulating 
regime, it is the localization length. Both diverge with the same exponent i^, 

From the scaling theory ^, Q , the conductivity exponent s and the length scale 
exponent v are related 

s = (d - 2)v (4) 

where d is the dimensionality of the system. This is called Wegner's scaling law. 
Numerical finite size scaling studies allow v to be accurately determined, while in the 
experiments s is measured. 

The most accurate way of estimating the exponent v is the transfer matrix method 
used in conjunction with finite size scaling||^, ^. For three dimensional (3D) systems, 
we consider a very long bar with cross section L x L. From the exponential decay 
length along the bar we define the dimensionless quantity 

Al = ^, (5) 
and assume a single parameter scaling form, 

Al = f{^/L) = F{SxL'/n- (6) 

where Sx = (x — Xc)/xc is the distance from the critical point. Assuming that A^ is 
an analytic function of x for finite L, we expand (^ as 

Ai = A, + aiLi/-(^ _ + a^lL^/^^ix - x,)]^ + ■■■. (7) 
Fitting the numerical data to (|^), we estimate A^., v, x^ and the expansion coefficients 

In some cases, such as quantum percolation where the lattice structure is dis- 
ordered, the transfer matrix method is difficult to apply and we have to use other 
quantities which obey the single parameter scaling law (^). The analysis of energy 
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level statistics is free from the lattice structure problem. For example, when we di- 
agonalize a L x L x L cube and define by integrating the nearest neighbor level 
spacing P{s) up to some point so 

Al - / " P{s)ds, (8) 
Jo 

then the same scaling law is valid. In the insulating regime, in the limit of large 
L, Al approaches Ppoisson(s)ds where Ppoisson(s) = e^'* is the Poisson distribu- 
tion. On the metallic side, it approaches PvvD('S)ds where -Pwd(s) is the nearest 
neighbor spacing for Wigner-Dyson statistics. By making full use of scaling behavior, 
the exponent v has been estimated for various models, and the universality of the 
transition confirmed with high precision. 



2.1 3D Anderson model 

The Anderson model has the Hamiltonian 

H= exp(i0,,,)C|Q + J2 W.Cja, (9) 

<i,j> i 

where C|(Ci) denotes a creation operator of an electron at the site i on the 3D cubic 
lattice and Wj denote the random scalar potential at the site i. If the transition is 
universal, even when the distribution function of the random potential is changed, the 
critical exponent should be invariant. To demonstrate this, we set all the phases 9ij 
to be zero so that the time reversal symmetry exists (orthogonal universality class), 
and consider three different types of random potential: the box distribution 

l/W i\W,\<W/2), 
( otherwise ) , 

the Gaussian distribution 

1 / 

with (T^ — W'^/12, and the Lloyd model in which Wi has a Lorentz distribution 
W 

For this distribution all moments higher than the mean are divergent and the param- 
eter W is proportional to the full width at half maximum of the distribution. 

We then analyze the case when time reversal symmetry is broken by magnetic 
fields. If the transition is universal, the critical exponent should be independent of 
how we break the time reversal symmetry. We apply uniform magnetic fields and 
random magnetic fields. In the former case, the phases of the hopping are chosen so 
that the fiux per unit square in the xy-plane is 1/3 and 1/4 of the flux quantum, while 
in the latter the phases of the hopping are random. In Table |l|, we summarize the 
results of the critical exponents for various cases of the Anderson model. From this 
table, we see that v = 1.57 ± 0.02 for the system with TRS. When the TRS is broken 
(unitary universality class) either by uniform or random magnetic fields, v becomes 
1.43 ± 0.04, clearly distinguished from that in the presence of TRS. 
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Ac 




ReL 


OG 
OL 

flux 1/3 
flux 1/4 
random flux 


21.29(27,32) 
4.27(25,28) 
18.316 ± .016 
18.376 ± .017 
18.80 ± .04 


yj.o t uyt o, i oj 

0.576(73,78) 
0.579(76,88) 
0.5683 ± .0013 
0.5662 ± .0016 
0.558 ± .003 


1.58(55,61) 
1.58(46,65) 
1.43 ± .04 
1.43 ± .06 
1.45 ± .09 
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Table 1 The best fit estimates of the critical disorder Wc, Ac and the critical exponent 
and their 95% confidence intervals. OB, OG, OL are for the system invariant under the 
operation of the time reversal with site potential distribution function box, Gaussian and 
Lorentz, respectively. In ref.|7|, corrections to scaling are taken into account. 



2.2 Other models 

If the Anderson transition is truely universal, the exponents should be the same, 
irrespective of whether we use the Anderson tight binding model or not. One example 
is the stack of quantum Hall layers, where electrons are allowed to hop from one layer 
to another |10, n|. The critical exponent v is estimated to be 1.45 ±0.15, ^ consistent 



with that in the case of Anderson model with broken time reversal symmetry. 

Another important example is that of quantum percolation where the Hamiltonian 

is 

H = Y.iU,CjC, +h.c), (10) 

(y) 

with the transfer integral 

, V exp{i9ij) ( for connected bond ), , . 



( for disconnected bond ) 



Here we consider the bond percolation problem. Bonds are randomly connected with 
probabilities p. 6ij is the Peierls phase due to magnetic fields. The underlying lattice is 
a three-dimensional cube of length L with periodic boundary conditions imposed. As 
we increase p, we reach a geometrical percolation threshold Pc and an infinite cluster 
is formed, however the wave functions are still localized and the system remains an 
insulator. If we further increase the probability p and reach pq, the wave functions 
become delocalized, and for p > pq current can flow through the system. This is the 
quantum percolation. Near the quantum percolation threshold Pq, the length scale 
diverges as ^ ^ \p — Pq\~'' . 

As mentioned above, this model is difficult to study with the transfer matrix 
method, so the energy level statistics are used instead [|l^. The estimate of is consis- 
tent with that of the Anderson model, 1.45± .11 with TRS (% = 0), and 1.25± .08 in 
the absence of TRS where 9ij is randomly distributed between — tt and tt. The critical 
exponent indicates that quantum percolation may be in the same universality class as 
the Anderson transition lfl3|. 



^We usually use one standard deviation for error bar except in Table |l[ 
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3 Scale invariance and boundary condition dependence at the transition 

3.1 Fractal dimension 

We now turn our attention to the properties just at the Anderson transition. It is 
well known that at the transition, the wave function shows multifractal structure 
[ p^ list [l7|, [T^, p^ , 20 1 which leads to the scale invariant behavior of conductance 
distributions[^|2g, § if H, || and the energy level statistics H, |^ |8[ |||, |^ |l| 



The direct way to investigate the wave functions is to diagonalize the Hamiltonian. 
This, however, is strongly constrained by the limited memory and CPU time. Instead, 
we calculate here the time evolution of wave packets to extract the information of 
fractal dimension [Q. We first prepare the initial wave packet |0) close to the critical 
point by diagonalizing a small cluster located at the center of the system. The time 
evolution of the state at time t is then obtained by 

\t + At) = U{At)\t) 

where U (At) is the time evolution operator. We approximate U{At) by a product of 
exponential operators 

UiAt) = e-'^^*/'* = U2{pAt)U2iil - 2p)At)U2{pAt) + 0(Ai5) (12) 

withp= (2-2i/3)-i and 

U2{At) = Q-^f^i^t/^^ . . . Q-iHg-iAt/2h^~iHqAt/n^-iHg-iAt/2h , , , ^-iHiAt/2h ^--j^g-j 

where Hi, • • • , Hq is a decomposition of the original Hamiltonian H — J^ , Hi in which 
each Hi is simple enough to be diagonalized analytically |39| . 
The square displacement of a wave packets is defined by 

r\t) = {t\r^\t). (14) 

In metallic phase, r'^{t) is proportional to Dt where D is the diffusion coefficient. In 
the insulating phase, it saturates to the square of localization length, At the critical 
point, the anomalous diffusion!^, ^ 

r2(i)^t2/<i^i2/3 (15) 
is expected. 

The fractal dimension D2 is estimated from the autocorrelation function 



C{t)^\ / dt'\{t'\0)\\ (16) 



t 

Since C(i) represents the inverse of the volume of the wave packet at time t, 

from (|l|) we obtain 

C{t) - t-^'/"^. (17) 
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^ 1 10 100 1000 t 10000 

Fig. 1 r^{t) vs. t at the critical point. Anomalous diffusion of ~ t^^^ is observed. 




Fig. 2 Autocorrelation function C(t) as a function of time. Due to the fractal structure of 
the wave function, the decay of C{t) is slow, ~ t~° *^, compared to t^^'^ expected for normal 
diffusion. 
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Eq. ( |l7| ) can be derived more rigorously in Fourier space pO| . 

In Figs. and ||, we show the resuhs of r^(i) and C{t) at the transition in the 
presence of a uniform magnetic field. The strength of the field is 1/4 magnetic flux 
per unit square lattice in the xy-plane. The Anderson transition takes place at Wc = 
18.41/. By diagonalizing a small cluster of 7 x 7 x 7 located at the the center of the 
system, we follow the time evolution of wave packets in 101 x 101 x 101 systems. The 
geometric average of C{t) over 10 potential configurations are performed. By fitting 
the data for t > AOh/V, the fractal dimensionality D2 is estimated to be 

D2 = 1.3 ± 0.2 

considerably smaller than the space dimension d = 3. The above value is consistent 
with the estimate in the case of a random magnetic field as well as layered systems 
in high perpendicular fields. P4I 



3.2 Universal distribution functions and its boundary condition dependence 

The above fractal structure at the transition leads to novel level spacing distribu- 
tion P{s) and conductance distribution P{g) at the critical point. Recently, it was 
pointed out that P{s) and P{g) depend on the boundary condition This was 

rather unexpected, since even completely different models show the same P{s) at the 
transition |Q. 

To investigate the origin of the boundary condition dependence, we study A^ with 
the fixed boundary condition (f.b.c.) imposed in the transverse direction of a long 
bar, and compare it with that in the case of periodic boundary condition (p.b.c). For 
simplicity, we concentrate here on the Anderson model with orthogonal symmetry. In 
Fig. H, we show the plot of A^ vs. the strength of randomness W for p.b.c. (+) and 
f.b.c. (o). The size L is 6,8,10 and 12, and the accuracy of raw data is 0.1%. The 
critical point should be indicated by a common crossing point in the data. We see 
that this is not clear for f.b.c. but seems to occur at « 15. 51^, while data for p.b.c. 
indicate W w 16.5V. This difference is physically unacceptable. In fact, a detailed 
study using ^-square fitting makes clear that the simple single parameter scaling (|^) 
fails. 

To overcome this, we have to introduce corrections to scaling |^, |^ to take account 
of surface effects 

Al = f{5wL^/'',b{W)L-') (18) 

with b{W) an analytic function of W around Wc and SW = {W — Wc)/Wc. We then 
expand this expression as 

Al = fbuiUSwL'/n + hidwL'/nKW)L-' + f2{SwL'/n^^^^^^^—^ + ■ ■ ■ -(19) 

Fitting the data to this expression neglecting higher order terms, we define the bulk 
part of Al 

Ai"'^ = huii.iSwL'^-') = KL~h{5wL^''')m)L-^-f2{5wL^'n ^^^^^!^ '^' -(20) 
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Fig. 3 Kl vs. W for p.b.c. (+) and f.b.c. (o). The size L is 6,8,10 and 12. 

Fig.^ shows the rcsuhs. Now a common crossing for f.b.c. is obtained at W/V — 
16.53 ± .22, and the critical exponent is estimated to be 1.61 ± .16, consistent with the 
p.b.c. case. 

The above results indicate that even in the limit of infinite system size, Aj?""* 
depends on the boundary condition. Aj?""* is 0.576 ± .001 for p.b.c. Q and 0.419 ± .03 
for f.b.c. This means that , when f.b.c. is adopted, the correlation length at the 
critical point for finite L is significantly smaller, leading to the fact that P(s) as well 
as P(g) tends to share more features of an insulator. In fact P{s) with f.b.c. is closer 
to a Poisson distribution |Q, and the mean conductance {g) is smaller when f.b.c. is 
imposed [|| ||. 



4 Summary and concluding remarks 

In this paper, we reviewed recent progress on the numerical study of the Anderson 
transition in three dimensional (3D) systems. Now the concept of the universality 
class is established numerically. 

In spite of the fact that we can perform simulations for larger sizes in two dimen- 
sions, the estimate of the exponent is not as accurate as in 3D except for the quantum 
Hall case |4^. For example, in the presence of spin-orbit scattering, the estimates of 
the exponent v in 2D are rather scattered e.g. 2.2 2.5 and 2.8 [|8[ ||], which 
can not be distinguished from the value of = 2.35 ± .03 obtained in the quantum 
Hall case. This may be due to the fact that the corrections to scaling is larger in 2D 
Further study is necessary to confirm the concept of universality in 2D systems. 
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Fig. 4 Same as Fig.^ but surface corrections are removed. 
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